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Abstract 

We present a systematic analytical approach to the trapping of a random walk by a finite 
density p of diffusing traps in arbitrary dimension d. We confirm the phenomenologically 
predicted e~ CdP * 1 time decay of the survival probability, and compute the dimension 
dependent constant Cd to leading order within an e = 2 — d expansion. 
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1 Introduction 



1.1 Motivations 

It has been over thirty years that the trapping of a random walk in a medium with 
absorbing traps is the focus of physicists' attention. The reason is twofold. First and 
more importantly the experimental relevance of modeling the diffusion of exciton 
has led to a compact formulation of the problem in terms of trapping of random 
walks. The second motivation is of theoretical nature. Much is known on the 
static traps case. Connections between trapping and Lifschitz tails in disordered 
electronic systems [jj]], statistics of rare events [^] and Brownian motion theory 
have been established. That the problem bridges to other areas of physics stems 
from the ubiquity of random walks and their applications. 

Here we will concentrate our study on the diffusion of a random walker in 
the presence of diffusing traps. The configuration of the traps evolves in time, 
instead of remaining frozen in its initial distribution. Despite the simplicity of 
its formulation, this consitutes a genuine many-body problem in which an infinite 
number of degrees of freedom are coupled: the positions of the traps relative to 
that of the random walker are correlated. 

Let us provide a more precise formulation of the problem we propose to inves- 
tigate. We are interested in the problem of a tagged walker (hereafter christened 
as "the walker") evolving in a medium in which a finite density p of diffusing 
traps is present. Both the traps and the tagged walker have the same diffusion 
constant. Here we have set the diffusion constant to | for aesthetic reasons. How- 
ever it would be of interest to investigate the diffusion constants ratio dependence, 
since, as shown in at sufficiently low trap mobility, there exists a rich crossover 
regime between static and mobile traps. When the walker and a trap meet on the 
same site then, between t and t + At, the walker has a probability (3dt to die. 



The static traps case was "solved" some twenty years ago using field-theoretic 
methods based on instanton calculus [Q, ||, ||] in the sense that the asymptotic 
behavior of the survival probability has been determined analytically. The one- 
dimensional case can be solved in the sense that the survival probability can be 
computed exactly by the standard tools of random walk theory [[7], 0J. In the mean- 
while the mobile traps case has resisted analytic approaches, even in one space 
dimension, and it is only very recently that extensive numerical studies were de- 
voted to unravelling its fine properties What can be inferred, however, from 
phenomenological arguments [11] or from low trap density expansions [JTJ], is that 
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the survival probability of the tagged walker decays with time S as 



Z(S) 



exp{-c dP S d / 2 ) d<2 
exp(-CdpS) d > 2 



(1) 



where Cd is a dimension dependent constant (universal in d < 2, but non univer- 
sal in d > 2). What makes the problem so analytically difficult to tackle? One 
answer is that we are dealing with a truly dynamical problem. There is no underly- 
ing static partition function, as opposed to the static traps situation. This prevents 
in particular, the use of quantum field theory methods which proved so successful 
for static traps. Of course other types of field-theoretic mappings exist (dynami- 



cal theories built after the Doi-Peliti mapping), but all of them fail [jlfj to predict 
the asymptotic behaviour of the survival probability: just as is the case for static 
traps, this is a strong coupling problem not accessible with the usual perturbative 
toolbox. Finally, it is interesting to note the close connection between the present 



trapping problem and the study of systems of vicious walkers [13] which exhibit 
some analytic similarities (the one trap problem is the two vicious walker problem). 

Our goal in this article is to cast the problem in a form suitable for a systematic 
analytic treatment. We will go beyond the existing results for the survival probabil- 
ity by setting our calculation in a systematic framework and by giving amplitudes 
to leading order in an e = 2 — d expansion. 



1.2 Notations 

Consider N random walkers in a volume V = L d and denote by p = y their aver- 
age density (eventually the thermodynamic limit N,V — > oo with fixed p shall be 
taken). The positions of those walkers are denoted by Xj(s), their initial positions 
Xj(0) being random in space. Let r(s) denote the position of the tagged walker, 
which starts from the origin at time 0. We employ Brownian motion functionals to 
describe the dynamics of the system. The action governing the dynamics of the set 
of N + 1 walkers is Aq + A[ nt , where 



A [r,{xi}] = - ds 



dr\ 



As) H« 



dxi 



ds 



(2) 



encodes the free motion of the particles. The interaction term reads 



A int [r,{x i }]=pJ2 / 5 d S ^)(r( S )-x 4 ( S )) 
„• Jo 



(3) 
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It describes the trapping of r by the mobile traps Xj's at a rate /3. 

We define two types of averages, that with respect to the free action with an 
index 0, 

<...) = J Prf«>(r(0))n / d -^-e- A ° W 

i 

and that with respect to the full interacting process, 

(...) = <...e-^')o (5) 
without any index. Normalization is chosen so that (l)o = 1, so that (1) < 1. 

1.3 Quantities of interest 

We define the survival probability after time S by 

Z(S) = (1) (6) 
and introduce the function Z(q, S) defined by 

Z(q,S)^(e^) (7) 

which allows to define the Fourier transform of the probability of presence a walk 
conditioned to survive up until time S: 

From the knowledge of G(q, S) it should eventually be possible to deduce the 
mean square displacement R 2 = —2d^ of a walk conditioned to survive. It is 
not a priori obvious whether anomalous diffusion should occur and investigating 
this issue is beyond the scope of the present work. 

1.4 Preliminary analysis 

1.4.1 Mobile traps, immobile tagged particle 

Here we give the exact result for the survival probability: denote by f(x, r) the 
probability that a random walk starting from x reaches the origin for the first time 
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exactly at time r, the survival probability of tagged particle (located at the origin) 
is given by 



Z(S) = V-» £ ^ f 1 - E /(x, r)) = (1 - 1 EE -)) 

iXi\i=l \ r=0 / xr 



(9) 



where N(S) (~ r (£)r( 2 -£) if <^ < 2 and ~ if d > 2) is the average of the number 
of distinct sites visited by a random walk after r steps (a functional derivation of 
this century old lattice result can be deduced from [14] : M(S) = f d d x((l — 

e -»/drfM (rW -« )))o = ^En>o ( " ir( ^%r r(e/2)r ) - This already re- 
produces the characteristic features of the conjectured behaviour. Note that the 
result @ is exact -a property already noticed in [|I5|]- as long as the tagged walker 
is immobile. 



1.4.2 Dimensionless variables 

It is interesting to note that d = 2 appears to be the upper critical dimension in this 
problem. We shall henceforth set e = 2 — d. This can be seen by scaling out the 
walk's length S. We build two independent dimensionless couplings 

u = 0pS, v = {A-K)- d / 2 p~ 1 S- d/2 (10) 

and further define 

w = uv = (4Tr)- d/2 (3S^ (11) 

The local trap density, 

N 

p(x,t)^E^ ) ( x - x ^ t )) w 
i=i 

is a Poissonian variable and it is therefore difficult to exactly integrate out its fluc- 
tuations. Note that diffusion noise is often well-described by a Gaussian white 
noise. Such an approximation allows density to become negative (if extremely low 
probability), and since those very regions of low trap density govern to dynamics, 
the importance of those unphysical configurations is enhanced which completely 
messes any analysis based on the Gaussian noise approximation. We are going to 
attempt an expansion in powers of the trapping probability (which is proportional 
to (3) and eventually take the limit in which this probability goes to 1. The trapping 
rate j3 can be viewed as an extra parameter that can be used to probe the model's 
properties. 
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2 Pertubation expansion for Z(S) 



2.1 General structure of the expansion 

We want to evaluate Z(q, S) = (e iq r ^^). First we expand e~ A,nt in powers of the 
trapping rate 0, which yields 



E 

n=0 



exp 



x exp 



x exp 



nl 



dsi...ds„ 



d d fci d d k n 



(2it) d "' (2ir) d 



JV 

E 

ii,...,i„=l 



/ dTy(r)X;kfe( Sj -r) 
-i V / dry i .(r).k J -e(s j - r) 

7~. JO 



(13) 



i=i 



We have denoted by y(s) = 4^ the velocity of the tagged walker and by = 
that of trap i Consider the n th order term. It involves an average over n 
random walkers x^,...^^. We depict a general n th order term by the follwoing 
diagram: In a general diagram, the number of horizontal lines stands for the number 
of intersections that have take place between times and S, this is also the power 
of (3 in the perturbation expansion. The number of vertical lines, at order n, stands 
for the number of distinct traps that will intersect the tagged walkers's trajectory. A 
diagram of order n involving m < n distinct traps will be proportional to u n v m = 
u m w n ~ m . It is not a trivial task to determine the proportionality constant; this 
will depend in a complicated manner on the "topology" of the diagram (the latter 
constant will be a function of e). 



2.2 An example: diagrams of order 4 

As an example, we explicitly compute the contributions arising from the term n = 
4 in the expansion (13). We depict those contributions by a series of diagrams: We 
now state the explicit evaluation from which one can induct the general properties 
of those diagrams, such as their leading UV divergence. Diagrams (a, b, c, d, g) 
exhibit no overlapping loops. 

(a) (14) 



6 



Xl 



x 2 



x 3 







o 

3 

CO 

era 



CD 



— 5 



Figure 1: The bold trajectory of the tagged walker r meets that of the traps xi, x 2 , 
X3. The depicted diagram is of order 6 (there are six intersection times). 
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Figure 2: Shown are diagrams of order 4. Diagram (a) is proportional to « 4 , dia- 
gram (b) to u 3 w, diagrams (c,d,e,f) to u 2 w 2 and diagram (g) to uw 3 . 



(&) =^ 3 -r(,/2) f ^ M (15) 



{c)= l -u\ w T(e/2)f^—^ (16) 



(d)= u 2 ( ^ r(e/2)) 2___ (17) 



Diagrams (e, /) possess an overlapping loop: 

(e) =u 2 w 2 I dxidx 2 dx 3 dx 4 ((x 4 - x 3 )(x 2 - x\) - \{%2 - x 3 ) 2 ) _ 2 

= y dt 1 dt 2 dt 3 6(l -h-t 2 - i 3 )(l -h-t 2 - t 3 )(*l*2 + *2*3 + «3*1 + 



1/50 824 , 4,2 - y/3\ . 

S -y- 481n3 + ^ ln2 "7! ln 2T71 



(18) 
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The coefficient | = ^jj^t, where D = \ is the trap diffusion constant. If the 
traps were static, one would have D = and one would recover a standard 4> 4 
propagator diagram with a double pole in 4 — d. A nonzero diffusion constant acts 
as a partial UV regulator which shifts the existence of a UV divergence from d = 4 
down to d = 2. 

(/) =u 2 w 2 ( dxidx 2 dx 3 dx4 (x 4 - 23)" (x 2 — X\ — j(x 4 - £ 3 ))~2 

/_d 3 d 

d£id£ 2 d£ 3 e(l - ti - t 2 - *3)(1 -h-t 2 - h)t 2 2 (h + -t 2 + t 3 )-2 

(19) 

Finally, 

The conclusion to be drawn from those explicit expression is the way the short- 
distance divergences (the e poles) are interwined with the topology of the diagrams 
which are considered. 



2.3 Nonoverlapping diagrams 

We now focus the subclass of nonoverlapping diagrams: in those diagrams, there 
are no overlapping loops. In order to fully characterize such a diagram, one needs: 

• its order, denoted by M; this is also the total number of intersections. 

• the number n of distinct traps intersecting the walker's path. 

• the number of times rrij trap j (j = 1, n) intersects (in a row) the walker's 
path. Of course, Y^j=i m j = M. 

• among the set {mj}" =1 , we count the number of m/s which are equal. If 
they take v distinct values, we call pt, £ = 1, v, the number of mj coef- 
ficients that take the same value indexed by I. This is necessary to properly 
determine the symmetry coefficient of each diagram. The allowed values of 
v run from 1 through v n = n if n can be written in the form n = k{k + l)/2 
(k an integer) and u n = n — 1 otherwise. 



9 



r xi 

I 



(2) 
'm2 



m n intersections 

/ 




n traps 



In principle the following calculation holds for mj > 2 but it can be seen to triv- 
ially extend to mj = 1 without modification. 

First consider a particular diagram in which the M intersection times ordered 

< < 4 1} < ... < sg\ < 4 2) < ... < S W < S (21) 

where the upper index j (j = 1, .., n) denotes the trap that intersects the walker's 
path and the lower index i counts the number of intersections of trap j (i = 
l...,rrij). Its graphical representation is as follows: We call Dm,u the value of 
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that diagram. Given that sequence of intersections, we have to evaluate 



Dm,h = 



Ml 

(exp 



1 



exp 



exp 



\- m - i).. ( jv - „ + 1)/ d.f ». .M& j ^r-^ 

/ dry(r) EE k ?W } -^) 
1/0 \j=l i=l / 

S f n m j 
J ° \ .7=1 i=l 



j=i i=i 



0) 



(_1)M _ !)_ - n + I) 
— Ml M™ 

rl rl—ti rl—ti...—tM-l 

x dti dt 2 ... / dt M 

Jo Jo Jo 

d d d_ ( 

J-0 j. 2 J- 2 J-O j. 2 4-0 j. I 

L 1 i 2 u 3 ... t mi+ i L mi +2 •■■ b mi+...+m n -i+l ••• b M 

{-l) M N(N - 1)..(N - n + 1) u n (wT(^)) M - n 
_ _ r(n + l + (M-n)| 



M! 



N n 



(22) 



While the explicit evaluation is not a trivial task, the final result is particularly sim- 
ple. 



We now consider a general diagram of order M with n intersecting traps char- 
acterized by {mj}j = i t ,,, tn , {pe}e=i,...,v We do not specify the order in which the 
traps intersect the walker's path. Such a diagram has a value 

Ml Jfr \ , „ M n! u»(wr$)) M - n 

W%^-WUVi\ n - [}}, j ) M ' n " ( j IILiP«'r(n + l + (M-n)§) 

(23) 

At a given order M, a diagram involving n distinct walkers has a higher e-divergence 
if it is free of overlapping loops. Retaining for each order M only the nonoverlap- 
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ping diagrams, we get the survival probability 



+00 M M 

^) = i+EE E E E (- 1 ) 

M=l n=l mi+...+m n =M v=\ p\+...+p v =n 



e > 
2' 



nLi^ !r ("+ i + ( M -^) 

(24) 

Since we are concerned with leading divergences, we will approximate 



1 



r(n+l + (Af-n)S 



1 



(25) 



Summing over all values of M — n yields a factor (1 + j^ry) n Setting g 
1+ Jr(£) ' we are ^ e ft w i tn 



= e~ 9 (26) 
In the limit u, w — > 00 with — fixed, this leads to the asymptotic behaviour 



Z(5) ~ e - 



(27) 



with c £ = 2ire + C(e 2 ). We emphasize that this is not a low density expansion 
since the variable — oc pS d l 2 is held constant. We could also have written 



Z(S) ~ exp 



'^L r(2 + n|) 

n=0 ^ . 



(28) 



which is, to leading order in e, equivalent to the previously mentioned result. 



The only approximation made is to neglect overlapping diagrams. This is 
equivalent to performing an e-expansion of c £ . Compare two diagrams that dif- 
fer by the presence of an overlapping loop instead of a non-overlapping loop: in 
the latter, the leading short distance singularity is not constrained by time ordering. 
In the former, the time ordering associated with an overlapping loop acts as a par- 
tial short-distance regulator. The related UV divergence will always be smoother 
than in the nonoverlapping counterpart. However the power of S that appears in 
Eq. ( j27| ) is exact since, at a given order in u and w, diagrams with or without over- 
lapping loops have the same S dependence. We argue that a diagram containing 
at least a pair of overlapping loops has a softer e divergence than the equivalent 
diagram in which the overlapping loops have been disentangled. 
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2.4 Perturbation expansion of Z(q, S) 

As before, we shall now keep in the perturbation expansion only diagrams contain- 
ing no loop overlap, this leads to 



Dm,u = 



(-/?) 



M 



Ml 



(exp 



cxp 



exp 



N(N-l)..(N-n + l) Jds?...ds%l J 



(2vr) d "' (2vr) d 



n rrij 



dry(r) [ qQ(S - r) + £ J>^6(^ - r) 

j=l i=l 



J ° \i=i i=i 



j=i i=i 



iV(iV-l)...(iV-n+l) M _ n 

"C — — (X U' 



M! " M n 

rl rl—ti rl-ti...-t M -l 

x dii / dt 2 --- / dt M 

JO JO JO 



i? t 2 2 e^ t 3 2 e^ ... tt^i tJ^>e ^~ ... t° 



9 2 s (-l) M 



=e 4 



Ml 



\ w T( £ -)) M ~ n J 



mi+l °mi+2 

dz 



d 1 St mi+2 
1 



2vri ( z + ^)n+i 



A q 2 St, 



mi+...+m„_i+l M 
(M-n)f 

(29) 



where the integration path runs from — ice + up to +ioo + 0. The final result now 
reads 



Z(q,S)=e- 3 T i [ ^ r°°d*e-^ + 4%x P 



ux 



1 + 



r( £ /2)W 2 



(30) 



so that G(q, S) = e 2 . Diffusion is not affected at this level of the approxima- 
tion. 



3 Discussion 
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3.1 Comparison with the static trap problem 

When traps are static, the Brownian motion functional depends on the tagged walk- 
ers trajectory only through its local time x(x, S) = J Q ds 5^(r(s) — x), hence 
it can be mapped onto an 0(n) field theory with an interaction term than can be 
shown to be 

A eff [r] = - J d d x In (l - p + ^/.'^W'M) -+-Jd d xhi(l-p + pe~P<f> 2 ) 

(31) 

We refer the reader to [|], ||, ^] and references therein for details on how to derive 
this correspondence and how to exploit it. For small values of (3 or of the field this 
produces an attractive — (p A theory (which has no minimum, and which is equiv- 
alent to truncating the walker's effective self-interation potential to that of a self- 
attracting walk —[3 2 p 2 J dsds' 5^(r(s) — r(s'))), hence the complete knowledge 
of the interaction potential is necessary. No truncation of the interaction potential 
can be used. 

Attempting an expansion in powers of (3 will generate self-interaction diagrams 
that must be evaluated in a standard way (a loop corresponds to a 5(r(s) — r(s')) 
interaction). Loops in diagrams occur when a trap is visited at least twice. All 
those features have their counterpart when traps are mobile. 

If traps are mobile, the best Langevin equation that describes the density of the 
diffusing traps is 

d t p-DAp = rj, (r]{x,t)r](x',t')} = 2DpVV'S(t-t')5 {d \x-x') (32) 

If this Langevin equation were exact (which it is not since the density is a Poisso- 
nian variable and not a Gaussian one), this would lead to an effective interaction 
term depending only upon the path r 

We now see that the 5^(r(s) — r(s')) -function appearing in the static traps case 
has been replaced by g(r(s)-r(s')) (with g(x,r) = (47rD|T|)- d / 2 e- x2 /( 4D l T D and 
D = 1/2). In fact, any given graph of the static traps problem has its mobile traps 
counterpart, which is evaluated simply by replacing the <5-function in the loops 
(whether overlapping or not) by the traps propagator g. There is a one-to-one 
correspondence between the graphs of the two problems. 



I\\2 



(r(«) ~ rQQ) 
4D\s 



(33) 
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3.2 Beyond the leading order 

_ „ d+1 

For static traps the mean-square displacement is known [|7|] to scale as R ~ S d + 2 

__ d 

and the return to the origin probability as S d + 2 . Within the framework of our ap- 
proximation, if traps are diffusing, we find that both the mean-square displacement 
and the return to the origin probability take their simple random walk asymptotics. 
Expanding the function G(q, S) beyond leading order in e shows that this picture 
might not hold for the full interacting theory. Indeed, the neglected UV divergen- 
cies could be responsible for a modified scaling behavior. For instance, one can 

show that G(q, S) = (l + ^F(q 2 S) + 0(1)), with F a known function 

that can easily be expanded in the vicinity of 0. Deeper insight into the structure of 
the perturbation expansion is required before exploiting such results in the spirit of 
the approach developped by Duplantier for polymers [|l6|]. 

3.3 Conclusion 

We have given the first systematic approximation scheme for the survival proba- 
bility of a random walk evolving in a medium infested with freely diffusing traps. 
Our method exploits the properties of low-dimensional Brownian motion. As op- 
posed to the static trap case, no mapping to an equilibrium-like field theory exists 
and due to the strong coupling nature of the problem the elegant field-theoretic 
methods of reaction-diffusion problems fail. Extending the Brownian motion for- 
malism to truly nonequilibrium processes therefore proved unexpectedly fruitful. 
The method which we have elaborated should now be put to the test on a series 
of related problems which we now list. The first question that comes to mind is 
concerned with the scaling properties of a surviving path: what is the scaling be- 
haviour of the walker's mean-square displacement between and S, provided it 
has survived up until that time. Questions of lesser importance concern the num- 
ber of distinct sites visited by the walker, its return to the origin probability or the 
algebraic area swept (still restraining the average to surviving paths). Numerical 
approaches are usually difficult to implement [^] when it comes to exploring those 
properties due to the poor statistics over surviving events. Future work should ad- 
dress those issues. 

Acknowledgments : The author would like to thank Henk Hilhorst and Bertrand 
Duplantier for repeated discussions on this work. 
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